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We consider the Pseudo-Gap (PG) state of high-T c superconductors in form of a composite order parameter 
fluctuating between 2p F -charge ordering and superconducting (SC) pairing. In the limit of linear dispersion 
and at the hotspots, both order parameters are related by a SU(2) symmetry and the eight-hotspot model of 
Efetov et al. [Nat. Phys. 9, 442 (2013)] is recovered. In the general case however, curvature terms of the 
dispersion will break this symmetry and the degeneracy between both states is lifted. Taking the full momentum 
dependence of the order parameter into account, we measure the strength of this SU(2) symmetry breaking over 
the full Brillouin zone. For realistic dispersion relations including curvature we find generically that the SU(2) 
symmetry breaking is small and robust to the fermiology and that the symmetric situation is restored in the 
large paramagnon mass and coupling limit. Comparing the level splitting for different materials we propose 
a scenario that could account for the competition between the PG and the SC states in the phase diagram of 
high-T c superconductors. 

PACS numbers: 74.40.Kb 74.20.-z 74.25.Dw 74.72.Kf 


I. INTRODUCTION 

Reflecting our rather poor understanding of the physics 
of cuprate superconductors, two kinds of theories are still 
debating whether the final solution for this problem will 
be a “bottom-up” approach based on a strong coupling 
theory 1-3 or rather a “top-down” approach, where symme¬ 
tries and proximity to a Quantum Critical Point (QCP) plays 
a dominant role 4-7 . The recently proposed Eight Hot Spots 
(EHS) model is a promising “top-down” approach to cuprate 
superconductors 8,9 . It reduces the Fermi surface to only eight 
points on the anti-ferromagnetic (AF) zone boundary and tak¬ 
ing long-range AF fluctuations between them into account. 
When the dispersion is linearized at the hot spots, one ob¬ 
serves surprisingly that an SU(2) symmetry relating the d- 
wave SC channel (Cooper pairing) to the d-wave bond or¬ 
der, or Quadrupolar Density Wave (QDW), (charge channel) 
is present. Moreover, an imposant pre-emptive instability (of 
order of 0.6 J, where J is the AF energy scale) in the form 
of a composite SU(2) order parameter emerges, that has been 
identified as a good candidate for the pseudo-gap (PG) state of 
those compounds 9-12 . Motivated by an impressive set of new 
experimental results 13-22 , this theory points out to the emerg¬ 
ing idea that charge order is most certainly a key player in 
the physics of cuprate superconductors, in addition to AF or¬ 
der, d-wave SC state and the Mott insulator phase. Angle- 
resolved photoemission spectroscopy (ARPES) experiments 
confirm as well the presence of modulations in the SC state 
for underdoped Bi2201 23,24 and has been interpreted either in 
terms of charge order or Pairing Density Wave (PDW) inside 
the PG phase 25-28 . Finally, we like to mention the recent in¬ 
terpretation of Raman resonances by a collective SU(2) mode 
which is the first experimental result that supports the idea of 
a composite PG 29 . 

The idea of an emerging SU(2) symmetry belongs to a wide 
class of theories which explain the PG phase of the cuprates 
through the notion of degenerate symmetry states between the 
d-wave SC order and another partner. Maybe the most famous 


attempt to such a unification has been the SO(5) theory which 
relates the d-wave SC state to the AF order 30 . Not less famous 
is the SU(2) symmetry which relates the SC d-wave order to 
the 7T-flux phase of orbital currents 1 . In both cases, the key 
question was to argue that the energy splitting between the 
two orders was small enough so that thermal effects would 
restore the symmetry above T c and below T*, which are the 
SC and the PG critical temperatures. The same question holds 
here for the SU(2) symmetry relating d-wave and QDW or¬ 
der. Although the EHS model in its linearized version ver¬ 
ifies the symmetry exactly, it is not clear if a more realistic 
Fermi surface, including curvature and the whole band struc¬ 
ture, will induce small or large energy splitting. Moreover, the 
EHS model relies on long-range AF fluctuations which medi¬ 
ate the interactions, but the experimental observation points 
out to short-range AF correlations which potentially will gap 
a whole part of the Fermi surface, as depicted in Fig. 1. A 
theory for “hot regions” instead of “hot spots” is thus needed. 

In this paper we address carefully all these issue by evaluat¬ 
ing the SU(2) splitting on realistic Fermi surfaces, for the two 
distinct components of the composite order parameter: the d- 
wave X -field in the charge sector (which forms the QDW order 
in the EHS model) and the d-wave A-field describing the SC 
pairing sector. We find that the splitting of the SU(2) symme¬ 
try increases with the mass of the paramagnons, but decreases 
with the strength of the coupling constant between AF fluctu¬ 
ations and conduction electrons. This opens a wide regime of 
parameters where the splitting is minimal - of the order of a 
few percents - and where the SU(2) symmetry is expected to 
be recovered through thermal effects in a regime of tempera¬ 
tures T c < T < T*. Above the PG temperature T* all traces 
of the short-range charge and SC field have disappeared. Of 
course the SU(2) symmetry holds in the whole temperature 
range between T c and T*, hence subleading charge instabil¬ 
ities which occur below T* do have their SU(2) partners in 
the form of Pairing Density Waves (PDW) 31,32 . We also study 
the effects of the Fermi surface shape in breaking the SU(2) 
symmetry - which is only preserved in the EHS model with a 


2 


7r 



perconductors in the first Brillouin zone of a square lattice. The order 
parameters spatially extent over hot regions, that are centered around 
the hotspots positions and which are coupled by the 2pp and the AFM 
coupling Q to opposed regions. 



QCP doping 


FIG. 2. Schematic phase diagram of hole-doped cuprate super¬ 
conductors. The effective mass of the paramagnon propagator m e ff 
serves as as measure for the distance to the quantum critical point 
QCP. In the SC phase m e ff is positive and vanishes at the QCP to be¬ 
come negative in the AFM phase. The bare mass mb ar e is defined far 
away from criticality at some higher temperature, indicated by the 
dashed line. 


linearized dispersion. We find that the splitting is small away 
from the points where X and A are maximal. In the physi¬ 
cal situation of a large paramagnon mass and also strong cou¬ 
pling, the maximum of X moves towards the zone edge lead¬ 
ing to bond order parallel to the x-y axes. All these findings 
point to the realization that, while being a secondary instabil¬ 
ity to AF ordering, charge order is a key player in the physics 
of the PG phase of the cuprates. 


II. MODEL 


We start from the spin-fermion model 8,9,31 with Lagrangian 
L = Ly + , where 


Ly — y/* (d?-\-£k?i (j) o) \f/ , (la) 

L 0 = i</.£»-> + ^((/. 2 ) 2 . (lb) 



FIG. 3. (Color online) Panel a): Variation of the maximum value 
of the gap functions \x\ and |A| as a function of the mass m. Note 
that in all materials the 2pp pairing in terms of \%\ vanishes abruptly, 
whereas the SC pairing in terms of |A| approaches zero asymptoti¬ 
cally when the paramagnon mass m is increased. Panel b): Variation 
of the level splitting as a function of the mass m. In both panels 
A =44. 


The fermionic field i \f describes the electrons which are cou¬ 
pled via to spin waves described by the bosonic field 0 . 
The effective spin-wave propagator is D~ l =y|co| + |q| 2 +m 
where m is the paramagnon mass which vanishes at the QCP 
and 7 a phenomenological coupling constant, which we esti¬ 
mate from its form in the EHS model 9 to be of the order 10 -5 . 
For notational reasons we also write q = (i(Q, q). Neglecting 
the spinwave interaction (u = 0 ) one can formally integrate 
out the bosonic degrees of freedom. In the spin boson model, 
this generates an effective spin-spin interaction of the form 


Sint — ^JqSgS- 


q 


( 2 ) 


It is convenient to use a fermionic representation of the spin 
operator S and consider in the following only the paramagnet- 
ically ordered phase in z-direction, so that 7=3/2/. The par¬ 
tition function then writes Z = f fZf'Fjexp (—So — Si nt ) with 


v> = InGoi'f',, (3a) 

k,G 

Sint = — ^ JqVk.o Vk+Q+qivVk'^dVk'-q-QiO- 

k,k' ,q,G 


(3b) 
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FIG. 4. (Color online) Panel a): Variation of the maximum value of 
the gap functions \%\ and |A| in YBCO as a function of the coupling 
A for different masses m. Panel b): Variation of the maximum value 
of the gap functions \% \ and |A| for different compounds as a function 
of the coupling A for the mass m = 0.5. 


where the bare propagator is 

Co k = diag(j£D £k, i (0 + £—k—p? ico £k+p?^H“£—k)? (4) 

and the spinor field = (V4,<r, Y-k- P ,s> Vk+ P ,a, V-k,s) T - 
Furthermore, J~ l = AD~ l /3 A 2 , a G {t,|} labels the spin, 

Q = (n, n) T is the AFM ordering vector and p stands for the 
2pp vector, as depicted in Fig. 1. Note that the chemical po¬ 
tential /I is implicitly subtracted from the dispersion £k. We 
select the SC and the 2pp channel by introducing the two order 
parameters 

A k = (vt a Y- k)S ), Xk = {Yk,aVk+p,a)- (5) 

The interaction S\ is now decoupled by means of a Hubbard- 
Stratonovich transformation. The partition function becomes 
(up to a normalization factor) 

Z = I ^MA, X }^p[-So-S hef f]. (6) 

The effective interaction is 




rUt 


%k%k+q~^Jq ^k^k+q 


k,q,G 


- l * i i&iPk, ( 7 ) 

k,G 


with k = k + Q and the matrix M is 


M k = 





—&k 

X-k 


( 8 ) 


The fermions in Eq. (6) can now be integrated out so that the 
partition function becomes 



with G 1 = G 0 1 — M. After functional differentiation of the 


free energy F = —T\nZ with respect to we obtain the MF 
equations in matrix form 

M, = Y^h-k'Gk'- ( 10 ) 

k' 


The matrix equation can now be projected onto the different 
components. We will consider here the case of two competing 
order parameters which can not be non-zero at the same point 
in k space. Therefore, we consider the equation for A with 
X = 0 and vice versa. The gap equations follow as 


A* = T £ J- k _, 


Aj, + e^, + c o' 2 


Xk = -XT £ J~ k _ 


Xk> 


G' ( ^ '® , - £ k ')- e k'+p) - Xy 


(Ha) 

(Hb) 


To solve these equations numerically, £k is parametrized 
in tight-binding approximation with the following parame¬ 
ters: YBCO 33 (parameter set tb2), Bi2201 34 , Bi2212 35 and 
Hg 1201 36 and for electron doped cuprates 37 . The momen¬ 
tum sums in Eq. (11) are then carried out by discretizing the 
k- space by rectangular and equidistant grids. To keep the nu¬ 
merical computations tractable we neglect the frequency de¬ 
pendence of x an d A 38 - The Matsubara sums are then car¬ 
ried out exactly in the limit T —* 0 and the momentum sums 
are performed over 200 x 200 points and over one Brillouin 
Zone (BZ). Moreover, note that the 2pp vector which con¬ 
nects two opposed FS points at ±pp depends on the external 
momentum k in Eq. (11) and is only properly defined on the 
FS. Since we expect that the main contribution to the momen¬ 
tum sum in Eq. (11) comes however from the hot-spot region, 
we make the approximation to take the 2pp vector constant 
and take the 2pp from the hotspot for arbitrary points in the 
first BZ, as depicted in Fig. 1. Throughout this article, we use 
10% hole filling (respectively 10% electron filling in the elec¬ 
tron doped case) and the bandgap is 10 4 K. To evaluate the 
strength of the SU(2) symmetry, we study the level splitting 
\X — A|/A. This parameter afford the study of the relative am¬ 
plitude between the QDW and SC order parameter, % and A. 
It vanishes for a perfect SU(2) symmetry and becomes closer 
to one for a complete SU(2) symmetry breaking. From the 
Fagragian we find that 0 A has dimension of an energy and 
(f) ~ m -1 / 2 . To estimate the coupling strength, we evaluate 
an effective energy via JE'eff = ~ m~ x / 2 X and away from 

criticality Eq ~ m 0 ' A where we take mo = 1 as bare mass, 
compare Fig. 2. 
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FIG. 5. (Color online) Gap functions \%\, |A| and \% — A|/|A max | (from up to down) for different materials in the first BZ. The compounds 
are (from left to right): YBCO, Bi2201, Bi2212 and Hgl201. The figures correspond to the small mass and small coupling limit (A = 40 and 
mass m = 10 -3 ; note that the bare coupling is small w.r.t. the renormalized one Eq ~ 40 K and E e ff = 1265 K) that is most unfavorable for the 
SU(2) symmetry. 



III. RESULTS AND DISCUSSION 

A. Mass and coupling dependence on level splitting 

In order to test the effect of the curvature on the level degen¬ 
eracy, we have plotted in Fig. 3a) the variation of the maxima 
of x and A with the paramagnon mass m for a fixed value of 
the coupling constant A. We observe a similarity between the 
various compounds that we have tested. In a wide range at 
low value of the mass, the SU(2) degeneracy between % and 
A is verified within a few percents. The existence of such a 
regime is an indication that a PG driven by SU(2) symme¬ 
try is possible in cuprate superconductors. As the mass is in¬ 
creased we progressively lose the level degeneracy with the 
parameter % abruptly dropping down while the paring A is 
asymptotically going down to zero when the mass increases. 
It is interesting to see that the SU(2) symmetry is weak for the 
electron doped and Fig 1201 compound which experimentally 
show much weaker signs of charge order 22,39 . We also find 
that the compound Bi2212 behaves sightly different than the 
other compounds in Fig. 3 and 4, although it is not clear at the 
current stage where this deviation comes from. In Fig. 3b) the 
level splitting is directly shown for all the compounds and the 
two regimes, the one at low mass where the SU(2) symme¬ 
try is obtained and the higher mass regime where \% — A|/A 


becomes of order one are clearly seen. Within the non linear 
cr-model associated to the present theory 9 , the SU(2) regime 
is a signature of the PG of the system, while the energy split¬ 
ting of the two levels is associated to the superconducting T c . 
We see that Fig. 3a) mimics the generic phase diagram of the 
cuprates where the PG line T* abruptly plunges inside the SC 
dome at some value of oxygen doping. 


Although it is very encouraging to see that the SU(2) regime 
has a non-zero probability to exist, one can wonder whether 
the paramagnon mass in hole doped cuprate superconductors 
is small, since typically the AF correlation length is of few 
lattice constants 40 . The issue is addressed in Fig. 4a), where 
the values of % and A are shown for a fixed mass as a func¬ 
tion of the coupling constant A. Here again a generic pattern 
emerges. For small A the SU(2) symmetry is broken, but sur¬ 
prisingly, above a certain threshold of A, the SU(2) symmetry 
is almost completely restored. As seen in Fig. 4a), the bigger 
the mass is, the stronger the coupling constant needs to be for 
the symmetry to be restored. Figure 4b) shows that this be¬ 
havior is quite general among the different compounds. Note 
however that the electron -doped compound is less sensitive 
to the effect of increasing the coupling constant, compared to 
the other hole-doped ones, for which the SU(2) symmetry is 
restored for large enough A. 
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B. Spacial dependence of the splitting 

The SU(2) symmetry is not only broken due to the curva¬ 
ture of the Fermi surface at the hot spots, but it is typically 
broken in the BZ away from the eight hot spots. Figure 5 
shows the typical shape of \%\ and |A| for four compounds 
under investigation for the most unfavorable case for the sym¬ 
metry, that is for small values of the mass and coupling con¬ 
stant. The level splitting is shown as a density plot in the 
bottom. It is rather small almost everywhere in the BZ and 
at the hotspot positions with maxima of the order of 20-40% 
around the “shadow” Fermi surface. The main learning from 
these plots is that the variations of the Fermi surface geometry 
gives a rather small departing from the SU(2)-degeneracy for 
a various range of compounds. In all cases, the SU(2) sym¬ 
metry is well respected at the hotspot positions. 

In Fig. 6 we place ourselves in the strong coupling and 
strong mass regime and plot the variation of \%\ , |A| and 
\X ~ A|/A. The level splitting is also shown in Fig. 6 and 
found to be much smaller than the previous case in Fig. 5, and 
has dropped to an order of 5-10% 41 . Interestingly, the typical 
shape of \%\ and |A| in the BZ has changed compared to Fig. 
5, with maxima now around the zone edge. This is the justi¬ 
fication that “hot regions” instead of “hot spots” is the correct 
description of hole doped cuprate superconductors within the 
spin-fermion model. Note that since the maximum of \x\ is 
now at the zone edge, the wave vector corresponding to the 
associated charge order is now parallel to the x/y axes of the 
system, in similarity with the findings of Ref. [12]. 


* A 



FIG. 6. (Color online) Generic picture of the gap functions \x\, 
|A| and \x — A|/|A max | in the first BZ for hole-doped cuprates, here 
explicitly shown for YBCO. The figures correspond to the large mass 
and large coupling limit (A = 160 and mass m = 0.5, so that Eq ~ 
160 K and E Q ff = 226 K) where the SU(2) symmetry is well respected. 


C. Global trends in cuprate superconductors 


The interplay between mass and coupling allows us to re¬ 
late global trends in the phase diagram for cuprate supercon¬ 
ductors with the strength of the SU(2) symmetry breaking. 
In the electron doped compounds the coupling between AF 
modes and conduction electrons is believed to be weaker than 
for the hole doped case. From the above analysis we find that 
the SU(2) symmetry is less respected in that case leading to a 
smaller PG dome and overall smaller SC, see Fig. 7. For the 
same reasons La compounds where the coupling is also be¬ 
lieved to be small behave similarly. On the other hand, hole- 
doped cuprates like YBCO live in the large mass and large 
coupling regime. This results in broad gapped regions in the 
BZ where the symmetry is well respected so that both the PG 
and SC dome are large, as shown in Fig. 6. 

Finally, let us mention that a major effect of a magnetic 
field is to invert the order of the level splitting between the 
CDW and the SC components 10,11 . This will favor the charge 
order compared to the SC pairing. We believe that Umklapp 
scattering can have a similar effect to invert the level ordering, 
but leave detailed investigation for further studies. 


FIG. 7. (Color online) Generic picture of the gap functions \x\, 
|A| and \% — A|/|A max | in the first BZ for electron doped materials. 
The figures correspond to the small mass and small coupling limit 
(A = 40 and mass m = 10 -3 , so that Eq ~ 40K and E e g- = 1265K) 
that is most unfavorable for the SU(2) symmetry. 


IV. CONCLUSION 

In conclusion, this paper gives firm ground to the intuition 
that the charge sector is a key player in the physics of cuprate 
superconductors. While the main instability is still the AF 
ordering, the d -wave bond order relates to the d-wave pair¬ 
ing through an SU(2) symmetry. We have shown that there 
exists a wide range of parameters where the SU(2) degen¬ 
eracy is fulfilled, which gives a natural explanation for the 
large PG regime observed in certain compounds. We argue 
that compounds like electron doped cuprates or the La com¬ 
pounds are outside the regime of SU(2) degeneracy, and the 
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more pronounced energy splitting is the reason for the weaker 
PG regime. 
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